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Abstract 

We present the solution of a linear Restricted Solid-on-Solid (RSOS) model confined to a slit. 
We include a field-like energy, which equivalently weights the area under the interface, and also 
include independent interaction terms with both walls. This model can also be mapped to a lattice 
polymer model of Motzkin paths in a slit interacting with both walls and including an osmotic 
pressure. 

This work generalises previous work on the RSOS model in the half-plane which has a solution 
that was shown recently to exhibit a novel mathematical structure involving basic hypergeometric 
functions 3^2- Because of the mathematical relationship between half-plane and slit this work 
hence effectively explores the underlying g-orthogonal polynomial structure to that solution. It 
also generalises two other recent works: one on Dyck paths weighted with an osmotic pressure in a 
slit and another concerning Motzkin paths without an osmotic pressure term in a slit. 

1 Introduction 

Solid-on-Solid (SOS) models describe the interface between low-temperature phases, originally in 
magnetic systems such as Ising-like models [TJ [2] , though now more generally. They are effectively 
directed models in d + 1 dimensions. The configurations involved in the linear (l+l)-dimensional 
case, modelling the interface in a two-dimensional system, have also been used to model the backbone 
of a polymer in solution [3j. The critical phenomena associated with this model describe wetting 
transitions of the interface with a wall [2]. For the SOS model the phase diagram contains a wetting 
transition at finite temperature T w for zero field and complete wetting occurs taking the limit H — > 
for T > T w [1]. The basic model is naturally described in the half-plane but it is also natural to 
describe it in the confined geometry of the slit. 

The linear SOS model with magnetic field and wall interaction was solved in [JJ. Recently the 
restricted SOS (RSOS), where the interface takes on a restricted subset of configurations, was solved 
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with the same interactions of field and single wall interaction in the half-plane [7] . This has proven to 
be mathematically quite interesting as both the method of solution and the functions involved were 
novel. It was found that the solution could be expressed as ratios of linear combinations of terms 
involving the basic hypergeometric function 3^2. Also recently the polymer models of Dyck paths [5] 
and Motzkin paths [8] in a slit with separate interactions with both surfaces have been considered, 
without field-like terms. Here the solutions in the slit prove interesting both mathematically and 
physically. They are of interest physically because the infinite slit limit was shown to be subtly 
different to the half plane, realising a separate phase diagram [5J. Mathematically the slit exposes the 
orthogonal polynomial structure of the problem and uncovers hidden combinatorial relationships [8]. 
Finally, Dyck paths in a slit with wall interactions and weighted by the area under the path, equivalent 
to a field term in the SOS models, have only recently been analysed [6j, and show a rich (/-orthogonal 
polynomial structure. To explore further this area of research here we consider the RSOS model in 
a slit geometry with both separate wall interactions and a field/osmotic pressure term in the energy. 
We derive the novel g-orthogonal polynomials for this problem which give us the exact solution of the 
generating function. 



2 The RSOS model in a slit 

The RSOS model we analyse can be described as follows. Consider a two-dimensional square lattice 
in a slit of width (or thickness) w > 1. For each column i of the surface a segment of the interface 
is placed on the horizontal link at height < r.; < w, and successive segments are joined by vertical 
segments to form a partially directed interface with no overhangs. The configurations are given the 
energy 

- PE = -K \n - n-! l-H^n + BoY^ <^,o + B w ^ S ri)W • (2.1) 

i=l i=l i=l i=l 

As in [7j, we discuss the RSOS model in terms of lattice paths. An RSOS path is a partially directed 
self-avoiding path with no steps into the negative x-direction and no successive vertical steps. To be 
precise, an RSOS path of length N with heights vq to has horizontal steps at heights n, . . . , rjy, 
and vertical steps between heights rj_i and r{ for i = 1,...,N, but no horizontal step associated 
with tq. This means that an RSOS path starts at height ro with either a horizontal step (if r\ = ro) 
or vertical step (if t\ 7^ ro), but must end at height with a horizontal step. Figure [1] shows an 
example. 

The partition function for the RSOS paths of length iV in a slit of width w > 1 with ends fixed at 
heights ro > and r^r > 0, respectively, is given by 



Zr(ro;ri] 



exp(-/3£ , (r ;ri)) , |r - r%\ < 1 
, |r - n\ > 1 , 



(2.2) 



and 



Z%(r ;r N )= exp(-fiE(r ;n, . . . ,r N )) , N = 2,3,..., (2.3) 

0<ri,...,rpf_i<w 
|rj-ri_i|<l 
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Figure 1: A typical RSOS configuration beginning on the surface and finishing on the surface with a 
horizontal step: each horizontal step is assigned a weight x, each vertical step a weight y, each unit 
of area a weight q, each step that touches the bottom surface an additional weight a while each step 
that touches the top surface an additional weight b. The length of the configuration shown here is 
N = 9 in a slit of width w = 3, and the heights are ro = 0, r\ = 1, r2 = r% = 2, = 3, r§ = 2, r§ = 1, 
rj = 0, r§ = 1, and rg = 0. The weight of this configuration equals x 9 y 8 q 12 a?b. 



where 

N N N N 

- f3E(r ;n,...,r N ) = -K^ \n - rvi| - H ^ n + B ^ d nt0 + B w ^ 6 n , w . (2.4) 

i=l i=l i=l i=l 

Here, we shall consider paths with both ends attached to the surface, i.e. we shall focus on the 
partition function 

Z^ = Z^(0;0). (2.5) 

We define 

y = exp(-K) , q = exp(-H) , a = exp(i?o) , and b = exp(B w ) , (2-6) 

so y is a temperature-like, q a magnetic field-like and a and b are binding energy-like variables, and 
write 

Z% = Z%(y,q,a,b) . (2.7) 

The (reduced) free energy is then 

k(w; y, q, a, b) = Jim — log Z%(y, q, a, b) . (2.8) 
Define the generalised (grand canonical) partition function, or simply generating function, as 

oo 

G w (x, y,q,a,b) = l+Yl xNz n(v, Q, a, b) . (2.9) 

N=l 

Thus, the radius of convergence x c (w; y, q, a, b) of G w (x, y, q, a, b) with respect to the series expansion 
in x can be identified as exp(— k(w; y, q, a, b)), hence 

k(w; y, q, a,b) = - log x c (w; y, q, a, b) . (2.10) 
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It is convenient to consider G as a combinatorial generating function for RSOS paths, where x, y, 
q, a and b are counting variables for appropriate properties of those paths. Interpreted in such a way, 
x and y are the weights of horizontal and vertical steps, respectively, q is the weight for each unit of 
area enclosed by the RSOS path and the x-axis, a is an additional weight for each step that touches 
the bottom surface while b is an additional weight for each step that touches the top surface. For 
example, the weight of the configuration in Figure [T] is x 9 y 8 q 12 a 2 b. 

If we send w — > oo then we recover the generating function of the half-plane 

oo 

G h P (x, y,q,a) = l+Y / x N Z^(y, q, a) (2.11) 
N=l 

where Z^f(y,q,a) = Z^(y,q,a,b) for w > N, noting that the paths can have no more vertical steps 
than horizontal steps in an RSOS path. 

We find easily the first few terms of G hp as a series expansion in x, 

G hp {x, y, q, a) = 1 + ax + {a 2 + ay 2 q)x 2 + ... , (2.12) 

where the constant term corresponds to a zero-step path starting and ending at height zero with 
weight one. 

Note that the radius of convergence of the half-plane generating function Xc P (y, q, a, b) is not a 
priori the limit of the slit x c (w;y,q,a,b) as was demonstrated in |5j for the corresponding Dyck path 
problem with q = 1. 

3 Exact solution for the generating function 

The key to the solution is a combinatorial decomposition of RSOS paths which leads to a functional 
equation for the generating function G w . This decomposition is done with respect to the left-most 
horizontal step touching the surface at height zero, and is shown diagrammatically in Figure [21 




G w (x,a) = l + axG w (x,a) + ax 2 y 2 qG w - 1 (qx,l)G w (x,a) 

Figure 2: The diagrammatic form of the functional equations for RSOS paths, indicating the combi- 
natorial decomposition of RSOS paths. 

We distinguish three cases: 

(a) The RSOS path has zero length, and there is no horizontal step at height zero. The contribution 
to the generating function is 1. 
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(b) The RSOS path starts with a horizontal step, which therefore is at height zero. The rest of this 
path is again a RSOS path. The contribution to the generating function is kxG w (x, y, q, k). 

(c) The RSOS path starts with a vertical step. Then there will be a left-most horizontal step at 
height zero, and removing this step cuts the path into two pieces. The left path starts with a 
vertical and horizontal step, followed by an RSOS path starting and ending at height one and 
not touching the surface, followed by a vertical step to height zero. This left path is effectively 
in a slit of width w — 1. The right path is again a RSOS path in the slit of width w. The 
contribution to the generating function is yxqG w -i(qx, y, q, l)ynxG w (x, y, q, k). 

Put together, this decomposition leads to a functional-recurrence equation for the generating function 
G w (x, y, q, a, b) = 1 + axG w (x, y, q, a, b) + aqx 2 y 2 G w -i(qx, y, q, 1, b)G w (x, y, q, a, b) , (3.1) 
with "initial condition" 

G (x,y, q ,a,b) = T -^-. (3.2) 
If we send w — > oo then we recover the functional equation of the half-plane 

G hp (x, y,q,a) = l + axG hp (x, y, q, a) + aqx 2 y 2 G hp (qx, y, q, l)G hp (x, y, q, a) . (3.3) 

whose solution is given in [7]. 
We can rewrite (|3.ip as 

G w (x, y, q, a, b) = - — 1 VT\ > ( 3 - 4 ) 

I - ax- aqx z y z G w -i(qx, y, q, 1, b) 

For general values of q, a formal iteration of (j3.4[) leads to a continued fraction expansion 

1 

aqx"y 



G w {x,y,q,a,b) = — . (3.5) 



1 — ax 



1 — qx 



q 3 x 2 y 2 



1 — q 2 x 



fx 2 y 2 



q 7 x 2 y 2 
T X 



bq 2w ~ 1 x 2 y 2 

1 - q w ~ l x 



1 - bq w x 

However, there is a non-trivial method to solve the functional equation for G as a ratio of q- 
orthogonal polynomials. It is clear that the generating function can be written as a rational function 

n i h s P w (x,y,q,a,b) 

G w {x, y, q, a, b) = —— — (3.6) 

Q w {x,y,q,a,b) 

though it does not simply follow to write expressions for these. It does however follow from the theories 
of continued fractions and orthogonal polynomials (see pages 256-257 of Andrews, Askey and Roy [9]) 
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that both the numerator P w and denominator Q w of the generating function satisfy recursions 



Pw(x,y,q,a,b) = < 



b, 

1 — bqx , 



w = 
w = 1 



(3.7) 



(1 - bq w x)P w -i(x,y,q,a,l) - bq 2w 1 x 2 y 2 P w - 2 (x,y,q,a,l) w>2 



and 



Q w (x,y,q,a,b) = < 



1 — abx , 

1 — bqx — ax(l — bqx(l — y 2 )) 



w = 
w = 1 



(3.8) 



bq w xQ w _ 1 (x,y,q,a, 1) - kg 2 " 1 1 x 2 y 2 Q w - 2 {x, y, g, a, 1) w > 2 . 



So it is the functions Q w (x,y,q,a,b) that are the g-orthogonal polynomials referred to in the Intro- 
duction. 

One can immediately note that 



so that 



P w (x, y, q, a, b) = Q w (x, y, q, 0, b) 
G w (x,y,q,a,b) ~ 



Q w (x,y,q,a,b) ' 

We now form the width generating function for the denominator as 



W(t, x, y, q, a,b) = J2 Q w (x, y, q, a, b)t w 

w=0 

and find a functional equation for W(t) from the recurrence (|3.8|) as 

W(t, x, y, q, a, b) = 1 — abx — abqx 2 y 2 t + tW(t, x, y, q, a, 1) 

— bqxtW(qt, x, y, q,a,l) — bq x y t W{q t, x, y, q, a, 1 



(3.9) 
(3.10) 

(3.11) 



(3.12) 



Unlike in the case of Dyck paths in a slit [6], one cannot solve for W(t) = W(t,x,y,q,a,l) by 
direct iteration, as the functional equation involves W(t), W(qt), and W(q 2 t). So let us return to 
Qw(x, y, q, a, b) and consider the case b = 1, as this is all that is required to find the full solution. Let 



R w = Q w (x,y,q,a, 1) . 
From the recurrence (|3,8p we have for w > 2 that 

flu, = (1 - q w x)R w - X - q 2w - 1 x 2 y 2 R w _ 2 

with initial conditions 

Ro = 1 — ax , 

R\ = 1 — qx — ax(l — bqx(l — y 2 )) . 



(3.13) 



(3.14) 



(3.15) 
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Attempting to mimic aspects of the half-plane solution [7J, we define S w via 

R w = (-l) w q w ( w + i y 2 x w S w . (3.16) 

This gives us the recurrence 

q w x(S w - S w ^ + y 2 S w . 2 ) + 5 w _i = (3.17) 

with initial conditions 

S Q = l-ax, (3.18) 
—qxS\ = 1 — qx — ax(l — bqx(l — y 2 )) . 

Continuing with inspiration from the half-plane solution [JJ, we try the Ansatz 

oo 

S w =fJ W J2 C n1~ nW ■ ( 3 - 19 ) 



n=0 



For n = we have 



and for n > we have 



This implies that satisfies 



xifj 2 - /i + y 2 )c = (3.20) 

(3.21) 



(y 2 q 2n — [iq n + fi 2 )qx 



fj 2 - [i + y 2 = (3.22) 
for our Ansatz to work. We now parametrise y via 

y 2 = A(l - A) (3.23) 

which changes (|3,20p into 

x(n-\)(iJt- 1 + A)cc = (3.24) 

and the recurrence into 



M Cn-x 



(3.25) 



(fi - Xq n )(fi - (1 - X)q n )qx 
We see immediately that either /x = Aor/i = l — A and that we have two solutions from our Ansatz. 
This leads to the general solution for S w as 

n-l 



oo / -. \ n n— i 

-n=n ^ / if— n ^ 



^)(A-d-A)^) 



(3.26) 



oo , 1 v n n-l u 

+B(1 - Ar Sb) ^'"Ri^wf- 



, , . , _ . A) - Xq k+1 ) ' 

n=0 x ' k=0 y * yvv 7 * y 

and so for via (|3.16p . One can then use the initial conditions (|3.15p to solve for the coefficients A 
and B. 
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Defining 



N 



pti ^(n— N —w+l)n 



i} :Q (<?; q)n(q; q)N- n (pq; ?)n(^; q)n- 



where 



n-l 



{t;q) n =\[{l-tq k ) 

k=0 

is the standard ^-product, after some lengthy calculations one finds 

Q w (x, y, q, a, 1) = <f (-+ 1 )/ 2 [(1 - ax)T^\x, q, A) + (1 - qx - ax(l - y 2 ))T^{x, q, A) 

with 

1 



A) 



1 - 2A 



^ q (N 2 -3N)/2^_ x ^ W -N 
N=0 



(1-A) 



A 



(1-A) 



X N- 



A 



and 

Ti, 2) (x,g,A) 



1 



1 - 2A 



7V=0 



? (iV 2 -JV-l)/2/_ ^w-l-iV 



w+1 V 1 — A ' 



(3.27) 



(3.28) 



(3.29) 



(3.30) 



(3.31) 



To obtain an expression for G w one can then substitute (|3.29p into (|3.8p to obtain an expression 
for Q w (x,y,q,a,b) and then this into f)3.10f) to give G w (x,y,q,a,b). 



4 The infinite width limit 



For any finite length the partition function in the slit becomes equal to the half-plane for large enough 
widths w. This means that the generating functions G w approach G p , when they converge. The 
half-plane solution can be found in [TJ. Here we take the limit w — > oo of the solution derived for finite 
width. This gives us a more compact expression for the denominator of the G hp than appears in [7J. 
After some work we find, using A = 1 — A, 

Q hp (x, y, q, a) = (1-ax) (x, q, A) + P (2) (x, q, A)] +(l-qx-ax(l-y 2 )) |>( 3 > (x, q, A) + P (4) (x, q, A) 

(4.1 



with 



oo oo 



pV( X ,q,\)=EE 



(— X ) M \ M + m A 1_m q^M A +Mm+rn 2 + ±M-rn 



M= 0m =o ( A ~ A)(g;g)oo(g;g)m(Ag/A;g) m (Ag/A;g)oo 

oo oo , \l-m \m+M 4\M 2 +Mm+rn 2 + \M-m 



(- X ) M A 1_m X r ' 



pM(x,q,\) = y: V 

M=0m=0 (A - \){q;q)oo{q;q) m (M/^q)oc{\q/\;q) 



oo oo 



P( 3 )(X,,,A)=EE 



{-x) M A M+m+1 \- m q^2 



M 1 +Mm+m 2 + ±M+m 



M-- 



om=o ( A ~ X){q\q)oa{q;q)m{M/^q)m{M/^q)z 



(4.2) 
(4.3) 
(4.4) 
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and 

00 00 ( rr\M \—m \M+m+l hM 2 +Mm+m 2 + hM—m 

PW (x, q, A) = V V ( 7 X) A A ^ ^ . (4.5) 

M=0m =o ( A - A)(g;g)oo(g;g)m(Ag/A;g) m (Ag/A;g)oo 

As in we then find 

^ftp/ \ Q hp (x,y,q,0) 

G p {x, y, q, a) = —r- . (4.6) 

Q hp {x, y,q, a) 

While this expression is equivalent to that found in [7J, we note that its structure is fundamentally 
different. 



5 Mapping to Motzkin paths polymer model 

As foreshadowed in the Introduction there is a mapping between RSOS configurations and Motzkin 
paths. Starting at the leftmost site one uses the mapping in Figure [3] to construct a Motzkin path 
from an RSOS configuration. For example, the mapping of configuration in Figure 1 can be seen in 



Figure 3: The mapping from RSOS paths to Motzkin paths. 
Figure [H This implies that the generating function we have found is also that of of Motzkin paths in 




/ \ 




Figure 4: The mapping of RSOS configuration in Figure 1 to a Motzkin path. 

a slit weighted by area (which can be used to model an osmotic pressure [6]) an interaction with both 
walls. In this way the current work is a generalisation of that found in [6j. 
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6 Conclusion 



In this paper we have presented a solution to the linear RSOS model in a slit with field and wall inter- 
action terms in the energy. In particular we have evaluated the generating function and demonstarted 
that its limit is the half-plane solution found earlier. The numerator and denominator polynomials of 
the slit generating function are novel (/-orthogonal polynomials associated with the continued fraction 
expansion of the half-plane solution. 
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